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INTRODUCTION 


llie  tropics  has  been  known  as  a  region  of  great  heat  for  many 
centuries.  The  Greek  division  of  the  earth  into  "klimata"  portrayed  a 
"torrid"  zone  between  the  Tropics  of  Cancer  and  Capricorn,  and  Hadley's 
eighteenth  century  (Hadley,  1735)  explanatic'n  of  the  causes  of  the  trade 
winds  was  based  on  the  idea  of  the  excess  heat  of  low  latitudes  causing 
equatorial  air  to  rise  and  produce  an  inflow  of  air  from  north  and  south. 
Among  more  recent  workers,  Siitpson  (Sinpson,  1930)  calculated  the 
atmospheric  heat  budget  and  showed  the  region  within  approximately  30* 
latitude  of  the  equator  to  be  one  of  positive  heat  balance  throughout 
the  year. 

Most  of  the  solar  radiation  which  is  the  source  of  all  this 
energy,  passes  through  the  atmosphere  to  the  surface  of  the  earth. 
Terrestrial  long-wave  radiation  causes  some  loss  of  surface  energy 
and  gain  of  heat  to  the  atmosphere,  but  the  greater  part  of  the  heat 
received  by  the  tropical  atmosphere,  and  passed  on  to  higher  latitudes, 

IS  contributed  by  the  turbulent  and  convective  processes  associated  with 
the  fluxes  of  latent  and  sensible  heat  from  the  tropical  seas  emd  lands. 
The  exchange  of  energy  at  the  tropical  air/surface  interface  is  thus 
extremely  important  to  the  global  atmospheric  circulation.  Mo-aover, 
these  transfers  are  frequently  small-scale  phenomena,  even  though  they 
ultimately  contribute  to  atmospheric  circulation  on  the  largest  scale 
(Kuettner  and  Holland,  1969).  Thus,  the  surface  exchanges  in  the  tropics 
need  to  be  studied  in  detail  and  on  a  microscale  if  their  contribution 


to  the  atmospheric  system  is  to  be  properly  understood. 

A  great  deal  of  attention  is  currently  being  paid  to  the 
air/sea  interface  within  the  tropics  (Gars tang  and  LaSeur,  1968;  Kuettner 
and  Holland,  1969) .  This  is  not  surpri.sing  in  view  of  the  Icurge  amount 
of  water  lying  within  30°  latitude  of  the  equator,  the  upper  layers  of 
vAiich  constitute  a  vast  storehouse  for  the  short-wave  solar  radiaticxi 
received  at  the  surface  (Kuettner  and  Holland,  1969) .  Hie  land  areas  of 
the  tropics  cannot,  however,  be  neglected.  Continental  territ'-ries  sucdi 
as  tropical  Africa,  tropical  Soutli  America,  and  northern  Austredia,  provide 
large  areas  of  distinctive  land-surface  energy  budgets,  vdiile  scattered 
through  the  tropical  oceans  are  numerous  islands.  The  latter  contribute 
to  the  atmospheric  heat  engine  in  two  ways:  by  the  energy  batlance  of 
their  :>urfaces;  and  by  the  turbulence  arising  from  their  nature  as 
obstructions  v’ithin  prevailing,  large-scale,  oceanic  wind-flow  systems. 

One  such  island  is  Barbados,  the  most  easterly  of  the  Leeward 
Islands  of  the  West  Indies,  13°  latitude  north  of  the  equator,  and  located 
well  within  the  trade  wind  belt  of  the  North  Atlantic  ocean.  There  is  no 
lend  between  Barbados  and  the  west  coast  of  Africa.  Thus,  the  island's 
166  square  miles  of  varied  surface  constitute  a  first  obstruction  in  the 
moist.  North  Atlantic  trade-wind  flow  after  many  miles  of  trans-oceanic 
passage.  Such  a  location,  in  itself,  provides  a  valuable  base  from  which 
f)  study  the  influence  of  a  heated  island  within  the  trade  wind  system 
(Garstang  and  LaSeur,  1968).  The  island  is,  however,  both  large  enough 
and  varied  enough  in  surface  character  to  provide  a  valuable  base  for 
sanpling  the  energy  budget  of  an  island  land  surface  under  tropical 
conditions.  In  the  south,  where  the  island  is  widest,  there  is  the 
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open,  rolling  country  of  the  St.  George  Valley  (Fig.  1) .  iBuediately 
to  the  north  the  land  rises  rapidly  and  relief  intensity  increases. 

Mont  Hillafay,  1,115  feet  aboVe  sea  level  and  the  island's  highest  point, 
is  rru^ly  the  centre  of  a  long  north-*south  ridge  %diich  divides  the  island 
Aazply  across  the  line  of  the  prevailing  trade  winds.  The  eastern  side 
of  th  ridge  is  ste^  and  deeply  dissected;  its  western  side  slopes 
saootiily  to  the  shores  of  the  parishes  of  St.  Janes  and  St.  Peter,  while 
gentle  slopes  also  diaracterise  the  terrain  descending  towards  Bridgetown. 
This  relatively  sisple,  yet  strong,  ^ysiografhy  thus  offers  sufficient 
variety  to  provide  a  base  for  stocking  how  relief  affects  surface  radiation 
and  energy  budget  characteristics.  At  the  sane  tine,  it  does  not  produce 
sc  coaplicated  a  pattern  as  to  develop  a  degree  of  diversity  idiich  defies 
analysis. 

A  nacageable  diversity  within  a  relatively  sisple  franevozk 
also  characterises  the  s'  .rfaoe  ei^ression  of  cultural  activities.  The 
cultivation  of  sugar  cane  dcuinates  the  agricultural  scene.  However, 
the  It^al  reqnizcnent  that  ten  per  cent  of  plantation  land  should  be  in 
food  crops,  ensures  the  existence  here  and  there  of  agricultural 
landscapes  other  than  sugar  cane  grcxring,  idiile  the  existence  of  Bridgetown,, 
together  with  scattered  areas  of  woodland  and  grassland,  offers  a 
total  enseid>le  of  cultural  landsc^es  sufficient  to  enehle  the  influence 
of  these,  as  well  as  relief,  to  be  incorporated  into  the  energy  budget 
investigations . 

This  ooirioinaticn  of  location,  relief,  and  cultural  influences 
thus  — Barbados  a  particularly  suitable  base  for  storing  the  surface 
energy  balance  of  a  tropical  island.  Such  a  stu^  involves  solving  the 
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I  basic  energy  balance  equations  which  for  a  land  surface  are: 


g 


=  (fi  +  q)  (1  -  a)  +  L+  -  Lf  (1) 

-  G  =  LE  +  H  (2) 

where  is  net  radiation,  Q  is  direct  short-wave  radiation,  q  is 
sky-diffuse  short-wave  radiation,  Lf  is  long-wave  emission  from  the 
atmosphere,  Lf  is  long-wave  emission  from  the  ground,  G  is  heat 
conduction  in  the  ground  surface,  and  I£  and  U  are  the  latent  and 
sensible  heat  fluxes  respectively. 

Both  these  equations  have  been  evaluated  by  numerous 
scientists  under  a  variety  of  conditions,  both  tropical  and  non-tropical. 
Die  instrumentation  and  techniques  for  their  solution  are,  therefore, 
well  known.  As  a  rule  the  work  has  concentrated  on  observations  at 
points  for  different  kinds  of  horizontal  surface.  Sudi  observations  have 
proxriuisd  an  insight  into  the  physical  processes  euid  diverse  influences 
involved  in  the  surface  energy  balance  and  earth/6d.r  interactions.  For 
a  realistic  picture  to  emerge,  however,  it  is  necessary  to  evaluate  also 
the  surface  variations  in  the  bcd.ance  of  energy.  Ohis  is  particularly 
i]:^rt:Uit  in  the  tropics  where  solar  energy  is  passed  into  the  atmosphere 
from  the  siurface,  and  where  smcdl-scale  exchanges  at  the  earth/air  inter¬ 
face  can  have  such  important  ccsisequences  for  the  circulation  of  the 
atmc-sphere  at  l£u:ge. 

It  is  connon  practice  in  climatology  to  evaluate  the  spatial 
distribution  of  phencmena  from  c^servations  taken  at  points.  Such  an 
evaluaticx:  is  as  accurate  as  the  basis  on  which  the  interpolaticxi  is  m'de. 


If  this  basis  is  wseJci  the  resiilting  patterns  of  surface  veuriation  are 
poor,  however,  accurate  the  original  point  observatians  mi^t  have  been, 
ihua  it  is  important  to  understand  the  processes  underlying  a  spatial 
diSbXibution  and  to  evaluate  the  surface  variations  in  these  terms 
from  point  observations,  rather  tlian  to  interpolate  blindly  from  point 
observaticAs  ttfhieh  may  only  ei^rcss  the  consequences  of  physical  processes 
and  which  are  not,  therefore,  in  theAiselves  a  valid  basis  for  drawing 
iseirithmic  lines  of  the  element  to  be  evaluated  in  spatial  terms  (Gemiier, 
1968) . 

It  follows  fztm  this  eurguomut  that  to  express  the  topographic 
variations  of  the  radiation  balance  (equation  (1} ) ,  or  cf  the  surface 
energy  budget  (equation  (2)),  it  is  necessary  to  see  these  as  egressions 
of  fundaanental  physical  processes  which  can  then  be  employed  in  spatial 
terms.  Both  equations  eipress  the  interplay  of  three  major  factors:  the 
influence  of  earth/sun  relationships;  the  effect  of  the  atmosphere,  largely 
expressed  throu^  weather  conditions;  and  the  contribution  of  surface 
(dxaracteristics.  Ihese  three  £acto3»  can  be  classed  respectively  cis 
macroscale,  mesoscale,  and  microscale,  in  both  a  spatial  and  a  temporal 
sense.  Thus,  earth/sun  relationships  auce  responsible  for  the  fundamental 
earthwcurd  flux  of  energy  which  varies  in  a  large-scale  way  with  latitude 
and  the  seauon;  atmospheric  influences  are  largely  mesoscale  in  chairacter 
in  terms  of  the  synoptic  scale  of  fundamental  weather  influences;  and 
surface  diaracter  operates  at  the  microscale  level  by  exercising  final 
control  over  conditions  at  a  given  spot,  within  the  framework  provided  by 
the  other  two  elements. 

E;q)eriroental  work  currently  preceding  in  Barbados  can  zu&aBuxe 
the  parameters  of  both  equaticsis  (1)  and  (2).  These  measurements  become 
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spatially  meaningful  when  they  are  applied  in  terns  ot  the  fundamental 

euproach  indicated  above.  This  provides  a  rational  way  to  solve  the 

radiation  balance  and  energy  budget  equations  in  terms  of  their  topographic 

Vcuriations.  Such  an  approach,  moreover,  does  not  necessarily  require 

a  close  network  of  observation  points  (Gamier,  1968) . 

Althea^  the  ultimate  aim  is  to  try  to  solve  both  equations 

(1)  and  (2)  in  this  way,  the  initial  concentration  has  been  on  the 

solution  of  equation  (1).  This  emresses  R  in  terms  of  shortwave  and 

n 

longwave  radiation  balances.  Ihe  work  to  date  has  eiqAiasised  the 
topographic  variation  of  shortwave  radiation  income  (Q  +  q)  the 
eveduation  of  tdiich  from  the  observations  of  a  single,  representative 
site  can  now  be  satisfactorily  achieved  (Gamier  and  Ohanira,  in  press). 

Ihis  evaluatic»  is  fundamental  to  the  solution  of  equation  (1)  and  also 
involves  the  most  intricate  calculations.  The  purpose  of  this  report 
is  principally  to  discuss  the  methods  used  and  to  show  some  of  the 
results  achieved.  Wbrk  on  the  evaluation  of  the  other  elements  of 
equation  (1)  is  currently  proceeding  and  will  be  reported  on  in  due 
course. 
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II 

THE  EVALOmON  OP  SOKFACB  VARIATIONS 
IK  SLObAL  RADIATION  INOOME 

In  order  to  evaluate  the  surface  variaticns  of  global  radiation 
incoae,  it  is  necessary  to  cud.calate  the  direct  and  s)cy-diffuse  ccnpunents 
separately,  and  tc  sun  the  results.  This  is  because  the  source  of  direct 
radiation  is  different  frost  that  of  sky-diffuse:  the  foraer  coses  directly 
from  the  sun  acting  eis  a  single  source  in  terse  of  its  position  in  the 
sky;  sky-diffuse  radiation,  on  the  other  heaid,  coses  fras  diverse  sources 
depending  partly  on  the  position  of  the  sun  in  the  sky  and  patrtly  on 
the  couposition  of  the  ataosphere,  especially  its  cloudiness,  at  a  given 
tlie.  Mhile  the  contribution  froa:  both  sources  can  be  treated  theoretically 
and  jHathenatically,  the  calculation  of  topogr^q^c  variations  of  direct 
radiation  is  slopler  than  that  of  sky^diffuse  radiation,  and  will  be 
considered  first. 

Evaluation  of  Direct  Solar  Radiation  Incone 


The  integration  over  all  wavelengths  of  the  fundamental  formula 
for  the  snnochroaatiu  intensity  of  the  solar  beam  yields  a  gener2u  law 
of  transmission  (Haltiner  and  Martin,  1957)  for  the  flux  of  direct 
solar  radiation  such  that: 

I^  -  IrP®  (3) 

where  I  is  the  energy  in  langleys  per  minute  delivered  by  the  sun  on 


a  surface  normal  to  the  sun's  rays,  is  the  value  of  extra-terrestrial 
radiation  equal  to  the  sol2ur  constant  divided  by  the  square  of  tiie  raditis 
vector  of  the  e^u^th,  p  is  the  mean-zenith-path  transmissivity  of  the 
atmosphere,  and  m  is  the  optical  air  mass. 

The  flux  of  direct  solcu:  energy  falling  on  a  surface  is, 
however,  cxily  equal  to  I  under  the  particular  situation  that  the  surface 

HI 

is,  in  fact,  normal  to  the  sun's  rays.  For  all  other  cases,  the  intensity 

of  rcidiation  falling  on  a  surface  is  the  value  of  I  modified  by  the 

m 

relation  between  the  angle  and  azimuth  of  the  surface  on  the  one  hand, 
and  the  hei^t  and  aizimuth  of  the  sun  on  the  other.  By  ejqtressing  these 
characteristics  of  the  surface  and  the  sun  by  means  of  unit  co-ordinate 
vectors,  an  equation  can  be  derived  to  express  the  flux  per  minute  of 
direct  radiation  on  any  surface  in  the  form: 

=  1^“  cosiXy^  S)  (4) 

where  I  is  the  intensity  of  the  shortwave  radiaticxi  in  langleys  per 
s 

minute  falling  on  the  sxirface,  X  is  a  imit  co-ordinate  vector  normal  to 
the  surface  and  pointing  away  from  the  ground,  ^  is  a  unit  co-ordinate 
vector  e3q)ressing  the  position  of  the  sun,  and  is  a  syntool  denoting 

the  angle  between  X  and 

This  equation  can  operate  in  several  ways.  The  way  chosen  in 
the  present  context  is  to  think  of  the  sun  as  moving  north  and  south  with 
the  seasons  along  the  meridian  of  solar  noon,  and  to  consider  the  surface 
as  rotating  around  the  earth's  polar  axis  once  in  twenty-four  hours.  This 
concept  enables  S  to  be  ejqjressed  solely  in  terns  of  the  cosine  and  sine 
of  the  sun's  angle  of  declination  (6),  and  X  to  be  expressed  in  terms  of 
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the  latitude  of  the  surface  (ft) ,  its  angle  of  slope  (given  by  the  zeni.th 
angle  (Z^)  of  30  ,  and  the  azinuth  (A)  of  the  surface.  It  can  then  be 
show)  (Gamier  and  CKiMura,  1968;  (Swura  1969)  that,  for  a  given  aOBent  of 
tine: 


cos  (jlyy 


Ksinil  cosH)  (-cosA  ::inZ^)  -  sinH(sinA  cosZ^) 


+  (cosjf  cosH)  cosZ^l  COSO  +  [cosi*  (cosA  sinZ^) 


-I-  sind  cosZ^]  sin6  (5) 

fdiere  H  is  the  hour  angle  neasured  frt»  solar  noon  positively  towards 
west  and  A  is  neasured  positively  fron;  north  throu^  east  and  negatively 
froB  north  throu^  west. 

Equation  (5)  indicates  that  the  value  of  cos(X^  £)  can  be 
derived  without  difficulty  front  surface  geometry  and  other  readily 
avcu.lable  data.  The  value  of  m  in  equaticm  (4)  can  also  be  derived  in 
terns  of  the  coi^tonents  used  to  e;q>ress  X  and  £  provided  that  the  secant 
j^roximation  for  the  value  of  m  can  be  used.  Under  these  conditions 
(Gamier  and  Ohnura,  1968) ; 

m  =*  1/cosZ  =  l/(cos6  cos|7  cosK  +  sin6  sin|7)  (6) 

s 

where  Z  is  the  zenith  angle  of  the  sun.  When  Z  exceeds  70®,  however, 
s  s 

the  secant  approximation  for  m  is  invalid  and  the  true  value  of  the 
optical  air  mass  as  given,  for  exanple,  in  the  Smithsonian  Meteorological 
Tables  (List,  1966,  p.  417)  must  be  used. 


-li¬ 


the  only  eleaent  of  eqnatlon  (4>  triiidt  camiOft  be  i^tainc^ 
fxoB  tables,  earth/sm  zelaticmshipe,  or  terraic  analyeia,  is  ataes^beric 
transadssivity  (p) .  this  is  a  highly  variable  qtunitit^,  beiag  in  a 
general  way  a  foncticn  of  the  weather  amd  prevailing  air  wass  characteristics, 
the  calcolation  of  p  for  cloadless  canditioas  can  be  made  by  standard 
■ethods  (de  Brlcfaaabant,  196.1).  However,  the  prooednre  for  evaluating 
direct  rzidiatian  suggested  here  assusees  the  exLstence  of  obistrvatiaiis  at 
a  site  representative  of  the  region  over  which  the  surface  variations  in 
short-wave  radiation  inooae  are  to  be  evaluated  (Gamier  and  Ohmra,  1968) . 


Ifeider  scdi  drciaBtanoes,  p  nay  be  evaluated  froa  equation  (3)  or  by 


using  the  eeasnred  value  of  direct  radiation  as  an  indicator  of  the  eean 
ataospheric  transaissivity  for  the  day  or  hour  in  question.  Both  these 
sethods  require  that  the  radiation  observations  being  used  pensit  direct 
radiation  to  be  distinguished  fron  sky-diffuse  radiation. 

Equation  (4)  is  tliie  fund£UBental  equation  for  cetlculating  the 
flux  per  minute  of  direct  radiation  on  a  surface  of  any  slope  and  azimuth 
anywhere  on  earth.  For  practical  purposes  the  equation  must  be  integrated 
over  time  periods  longer  than  the  minute  to  <rtiich  it  refers.  Ihis  integration 
can  be  achieved  by  means  of  a  summation  approximation  (Ohmura,  1968) 
using  sufficiently  small  time  intervals  (t)  in  the  expression: 


^st  "  ^r  ^  ^ 


where  is  the  total  over  a  given  tirje  period  such  as  a  day  or  an 
hour,  and  t^^  and  are  the  first  and  last  times,  expressed  in  hour 
angles,  when  the  sun  shines  on  the  surface  during  the  period  in  question. 
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Suinnation  at  one-minute  intervals  is  practicable  with  the  use  of  hi^- 
speed  conputers.  However,  such  a  refinement  is  not  normally  necessary. 
Analysis  shows  that  to  use  i.  =  20  minutes  will  give  ein  accuracy,  at  least 
for  daily  totals,  within  5%  of  that  using  t  =  1  minute  ((^imura,  1969). 
Since  an  accuracy  of  5%  is  within  the  observational  error  of  instmiments 
normally  used  for  measuring  radiation,  such  a  time  interval  has  been 
adopted  in  the  calculations  made  in  the  present  context. 


Evaluation  of  Sky-Diffuse  Radiation  Income 


The  evaluation  of  soleur  radiation  income  is  a  relatively 
stredghtforward  matter  in  that  the  soiurce  of  origin  of  the  radiation  is 
a  single  point  in  the  sky.  When  sky^diffuse  radiation  is  considered,  however, 
the  matter  becomes  more  con^licated.  Sky-diffuse  radiation  reaches  the 
surface  by  scattering  and  by  reflection  from  clouds.  Under  completely 
clear  sky  conditious.,  sky-diffuse  radiation  reaches  the  surface  only  by 
scattering,  with  a  maximum  concentration  of  the  source  of  diffuse  radiation 
in  the  area  around  the  sun.  The  presence  of  clouds,  however,  introduces 
additional  sources  of  sky-diffuse  radiation  xmtil,  under  overcast  conditions, 
the  source  of  sky-diffuse  radiation  is  a  homogeneous  celestial  hemisphere. 
Under  the  latter  conditions  the  scattering  by  water  droplets  will  be 
independent  of  both  wavelength  and  direction;  under  completely  clear 
conditions,  however, when  Rayleigh  scattering  dominates,  the  flux  of 
diffuse  radiation  depends  both  on  wavelength  and  the  direction  of  the  ray. 

It  is  clear  from  the  foregoing  that  an  almost,  infinite  number 
of  situations  exist  under  which  sky-diffuse  radiation  reaches  the  surface. 
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To  solve  all  the  cases  is  an  almost  impossible  task,  although  the 
mathematical  argument  involved  is  not  too  coo^licated.  Diis  argument 
may  be  e^qpladned  by  reference  to  Fi.q^  2. 

If  the  intensity  of  diffuse  radiation  from  the  part  of  the 
sky  (g,  H  )  is  expressed  as  D(g,  H  ) ,  the  amount  of  diffuse 
radiation  emitted  from  the  area  cosg  dg  d  H  and  received  on  a  surface, 
the  plane  of  which  is  norm2d  to  the  source  of  emission,  will  be: 

D(g,  )cosg  dg  d 

The  flux  of  diffuse  radiation  on  a  slope  of  gradient  0  will  then  be; 

D(g,  @  )cos^g  cos(  0+  0)dg  d  @ 

By  integration  it  can  then  be  shown  (Ohmura,  1969)  that  the  flux  of 
diffuse  radiation  from  the  vdiole  celestial  hemisphere  will  be: 

ir/2-e  Tr/2  ^ 

D  ■  /  /  D(g,  (li)  )cos“g  cos  (  (h)  +  0)dg  d  (h)  (8) 

®  -tt/2  -■ir/2  ^  W  ^ 

Equation  (8)  is  thus  the  fundamental  equation  for  calculating  the  flux 
of  sky-diffuse  radiation  on  a  given  surface.  To  solve  it  requires  a 
leurge  number  of  measurements  since  the  value  of  D(g,  0  }  changes 
rapidly.  However,  if  a  homogeneous  celestial  hemisphere  is  assumed, 
equation  (8)  reduces  (Ohmura,  1969)  to  the  well-known  Kondrat'yev 
equation  (Kondrat'yev,  1965)  for  calculating  the  sky-diffuse  radiation 
falling  on  a  slope  from  the  value  measured  on  a  horizontal  surface: 


where  is  the  sky-diffuse  radiatico  measured  on  a  horizontal  surface. 

This  eaxiation  is  simple  to  operate  and  uses  readily  obtained 
data.  For  this  reason  alone  it  tends  to  be  coaaonly  used.  In  the  ptesemt 
context,  moreover,  its  use  <^oes  not  lead  to  serious  inaccuracies: 
firstly,  if  the  formula  is  used  for  daily  totals  of  dcy-diffuse  radiation 
it  means  that  there  will  have  been  integration  over  the  entire  keni^pheze 
above  the  plane  of  the  surface,  so  that  tlie  inaccuracies  inherent  in  tiie 
formula  will  be  reduced  (Sobinson,  1966) ;  secondly,  the  contribution  of 
sky-diffuse  radiation  to  global  radiation  is  minimal  isder  clear-sky 
conditions  when  the  assumption  of  a  hanogeneoas  celestial  hemisphere  is 
least  VcJ.id,  and  is  maximad  under  overcast  or  very  cloudy  conditions  when 
calculation  by  the  formula  contributes  itu  moet.  accurate  results. 
Consequently,  in  the  present  study,  the  topographic  vanations  is  global 
radiation  are  calculated  by  combining  equations  (4)  and  (9) . 

Accuracy  of  the  Method  of  Evaluation 

Theoretical  consideratioiB  suggest  that  the  procedures  described 

for  evaluating  the  topographic  variations  of  global  radiation  income 

by  the  use  of  equations  (4)  and  (9)  will  yield  results  which,  at  least 

for  daily  totals,  do  not  exceed  the  5%  error  suggested  as  permissible  in 

view  of  instrumental  Invitations.  To  test  this  claim  a  series  of 

measurements  was  made  at  Mont  St.  Hilaire,  Quebec,  (lat.  45 "N)  from  August 

* 

29  to  September  2,  1968  . 


A  more  extensive  series  of  measurements  is  at  present  currently  in  progress 
in  Barbados,  but  the  results  from  them  are  not  yet  avfdlable. 


nzee  Kipp  and  Zcoen  p^anoaetexs  were  installad  at  an  angla 
of  20*  facing  north,  east  and  soittii  xe^ctctxw.lj.  At  the  sane  tine  an 
Bpplsy  pfranoneter,  tilted  at  an  angle  of  20*  on  a  tripod  and  suitably 
j^aded,  was  rotated  sucoessisely  every  twenty  Baamtes  to  face  north, 
east,  and  south,  iwo  other  Kipp  and  Zcoen  pyraanneters  recorded  global 
and  sky-diffase  radiation  ocntinoosly  on  a  horisontal  surface.  Ihe 
value  of  direct  shortwave  radiation  thus  evaluated  for  a  horizontal 
surface  was  used  to  obtain  the  transnission  coefficient  of  the  atnoephere 
uhich  was  tiien  eoployed  in  calculating  global  radiation  on  the  sloping 
surfaces.  For  this  purpose  the  following  equation  (Ohnura,  1969)  was 
used: 

Oo^. Global  -  I^**coe(Xy^  S)  +  E^cos^  j  +  A1  (Q  -  E^)8in^  ^  (10) 

In  this  equation  ,  the  first  and  second  terns  cm  the  righthand  side  are 
frcMB  equations '  (4)  and  (9) ,  Hhe  third  tern  is  reflection  fron  the 
surrounding  aurea.  Calculation  shows  that  inis  third  tern  is  generally 
less  than  3%  of  the  shortwave  reflected  redi.ation  and,  therefore, 
<K>ntribated  less  than  1%  of  the  global  radiation,  since  the  erperinent 
was  cmnductad  in  grassy  surroundings  where  the  albedo  was  less  than 
0.30.  Thus,  althou^  the  term  was  considered  during  the  experiment  its 
small  contribution  can  be  normally  safely  neglected  in  evaluating  surface 
variations  of  short-wave  radiation. 

The  results  of  the  study  are  shown  in  Fig.  3.  The  dots  are 
plotted  in  respect  of  the  twenty  minute  totals  conputed  and  mecisured 
during  the  four  days  of  the  ej^riment.  The  1j1  relationship  line 
on  each  diagram  indicates  visually  the  resulting  close  correlation. 


Fig.  4  Ccaputed  and  Neaaxired  Values  o£  Sky-Diffme 


Radiation  Values,  Mont  St.  Hilaire,  Quebec 


(lat.  45*N),  Aug.  29  -  Sept.  2,  l^)6& 
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statistically  this  rorrelaticn  is  calculated  as.; 


North  20® 
East  20“ 
South  20“ 


r  =  0.995 
r  =  0.996 
r  =  0.996 


These  results  suggest  that  the  aethod  for  calculating  topographic 
variations  in  short  vave  radiation  incone  suggested  in  this  report  is 
reasonably  accurate.  It  is  interesting  to  note,  however,  that  the  error 
in  calculated  results  is  greatest  on  both  north  and  south- facing  slopes 
%dten  measured  shortwave  radiation  values  axe  hi^,  but  that  the  sign 
of  the  error  on  the  north  slope  is  opposite  to  that  ontiie  south  sl<^. 

Ihis  can  be  explained  mainly  by  the  formula  used  to  evaluate  sky-diffuse 
radiati(»i.  A  comparison  of  the  computed  and  measured  vailues  of  this 
element  alone  is  shown  in  Fig.  4.  This  coopcurison  indicates  that  the  two 
vcdues  were  close  to  each  other  vdien  the  sky  was  overcast,  and  diverged 
considerably  vmder  clear  sky  conditions.  Thus,  when  the  sky  is  clear 
and  the  sun  is  hic^,  the  sky-diffuse  radiation  will  be  underestimated  cm 
a  slope  towards  the  sun  owing  to  the  high  proportion  of  sky-diffuse 
radiation  coming  from  the  vicinity  of  the  sxin.  For  the  same  reason, 
the  sky-diffuse  radiation  will  be  overestimated  on  a  slope  facing  away 
from  the  sun.  On  the  other  hand,  the  more  nearly  hoomgeneous  conditions 
of  low  sun  or  overcast  conditions  is  expressed  in  the  greater  accuracy 
of  the  computed  values  when  global  radiation  values  are  low. 


F=I.7J5« 
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Fig.  5  %e  Flow  Diagram  £or  the  Basic 
Coaipv^er'‘1Progra]nnie  to  calculate 
Direct  Radiation  on  a  Slope 


CX)14PaTER  PROGRAMMING  AMD  MAPPING  PROCEOORES 


'Ole  computer  programning  and  mapping  procedures  idiereby 
evaluati<»i  of  surface  variations  of  global  radiation  are  effected 
be  considered  in  Hiree  sections:  the  basic  computer  programning 
involved;  the  way  to  map  the  data;  and  some  particulau:  aspects 
useful  in  a  tropicad  context. 


Basic  CcMpputer  Progrannning 


The  conputer  programme  which  is  fundamental  to  all  calculation 

e 

and  mapping  of  global  radiation  on  slopes,  is  the  progrcunne  for  confuting 
direct  radiation  from  equation  (4) ,  using  the  elements  shown  in 
equations  (5)  and  (6)  respectively.  Hiis  progranrae,  arranged  for  a 
print  out  in  table  form  (Table  One)  is  illustrated  in  Figs.  5  and  6. 

The  notation  used  in  the  programme  is  as  follows: 

F  is  ths  latitude,  given  in  radians  . 

FO  is  the  emgle  of  declination  of  the  sun,  also  given 
in  radiains. 

SOL  is  the  value  of  the  extra-terrestrial  radiation  in 
langleys  per  minute. 

G  is  the  atmospheric  transmissivity,  given  in  the 

programme  illustrated  for  numbers  between  0.40  and 
0.85  in  steps  of  0.05. 

Q1  is  the  value  of  the  optical  air  mass;  this  is 

calculated  from  equation  (6)  imtil  the  sun's  zenith 
angle  is  equal  to  70®,  from  which  point  onwards  the 
true  optical  air  mass  as  given  in  the  Smithsonian 
Meteorological  Tables  (List,  1966)  is  used. 
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0001 

0002 

0003 

0004 

0005 

0006 

00C7 


OOCB 

0009 

0010 

0011 

0012 

0013 

0014 

0015 

0016 

0017 

0018 

0019 

0020 

0021 

0022 

0023 

0024 

0025 

0026 

0027 

0028 

0C29 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 

0038 

0039 

0040 

0041 

0042 

0043 

0044 

0045 

0046 

0047 

0048 

0049 

0050 

0051 

0052 

0053 

0054 

0055 

0056 

0057 

0058 

0059 


C  CALCULATION  OF  OIAECT  INSOLATICN  ON  THE  SLOPES 
C  UNIT  LANGLEYS  PER  OAV 

CIPEKSICA  X(3ti  01191 
C  LATITUDE  IN  RADIANS 
F-0.22t9 

C  CECLINATION  OF  THE  SUN 

130  READI5tl31. END-999 iFC.  SCL 

131  FORNAT  (F5.0tF4.0) 

C  TRANSNISSIVITV 

8  G-0*40 

7  NRITE  (6.101  F.  FG.  G 

10  FORNAT  (lH1.9HLATITUDEi.FT.4//13H  DECLINATION*. F6.2// 

116F  TRANSNISSIVITV*. F5.2/// 

211H  ANGLE  CF  .4SX. 

3  18HAZINUTH  OF  SLOPE  /6H  SLOPE. 3X.3H  0.3X.3H  10.3X.3 

4H  20.3X.3H  30.3X.3H  40.3X.3H  50.3X.3H  60. 3X. 

53H  70.3X.3H  P0.3X.3H  90.3X.3H1C0.3X.3H110.3X.3H120.3X. 
63H130.3X.3H140.3X.3H150.3X.3H160.3X.3H17C.3X.3H180/I 
8-0.0 
BA-0.0 
2  A-C.O 

DO  1  1*1.19 
Xdl— CCS(A|4SIN(BI 
XI2)*SIN(B|4SIN(AI 
X(3)*CCS(B) 

T1*(X(1J*SIN(F»*X(3»*CCSIFI|4C0S(F0) 
T2*(-X(ll*C0S(FI*X(3l*SIN('rl)*StN(F0t 
W— 2.1C 
01  11*0.0 

5  QC*C0$lFCI4CCS(FI*C0Stkl«SIN<FCI*SIN(F( 
h*h«0.C873 

IF(OO.LE.O.Oi  GC  TO  5 

02— X(2J*SIN(HI*CCS(F0»*T1*C0S(W» 

CT*02*T2 

IF(I<.GT.2.1I  GC  TC  1 
IF(QT.LE.O.Ot  GO  TO  5 
C1«1.0/(CC«0. 00011 
IF(Ql-2.8)  121.  121.  122 
122  IF(Ql.GE.2.9i  01-2.9C 
IF(01.Gg.3.0)  Gl«3.05 
IF(Q1.CE.3.1I  01*3.21 
IF(01.GE.3.3I  01*3.39 
IF(gi.C-E.3.5)  01*3.59 
IF(01.GE.3.7I  01*3.82 
IF(01.CE.4.0)  01*4.07 
IF(01.CE.4.6I  Cl*4.72 
IF(01.CE.5.0)  01*5.12 
IF(01.GE.5.5)  01-5.CC 
IF(Q1.GE.6.0)  01-6.18 

IF(01.GE.  6,5)  01*6.88 
IF(01.GE.  8.C!  01-7.77 
IFCOl.GE.  8.81  01-8.90 
IF(Ol.CE.lO.C)  01-1C.39 
IFIQl.GF. 13.01  01-12.44 
IF(01.GE.17.0)  01-15.36 
IF(01.GE.24.C)  01-19.79 
121  IF(ABS(C1I.GT.45.0)  C1-30.C 
Q(II»0(I)+20.04SCL-G*4CIA0T 
IF(H.LT.2.1)  GC  TO  5 
1  A-A*0. 17453 

kRITE(6,20l  BA,  0 
20  FORMAT (20F6.1) 

B-B-fO.  17453 
BA-BA+IO.O 

IF(6.LT.1.60I  GO  TC  2 
G-G«0.C5 

IF(G,LE.0.85)  GC  TO  7 
GO  TO  13C 
999  STCP 

fnd 


Fig.  6  The  Computer  Programme  to  Calculate 
Direct  Radiation  on  a  Slope 
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Table  One 

Table  of  Direct  Radiation  on  Slopes,  Barbados 


I 


L4TITU0E1  0.22(9 
0ECL1((TICN*  C.4C 
TR«NSMISStVITY«  0.6S 


ANGLE 

SLCPE 

OF 

0 

10 

20 

30 

50 

50 

60 

70 

AZHUTH  CE 
80  90 

SIC8E 

100 

110 

120 

l30 

150 

130 

160 

170 

ISO 

0.0 

512.; 

512.5 

512.5 

512.5 

512.5 

512.5 

512.5 

512.5 

512.5 

5l2.» 

512.5 

512.5 

512.5 

512.5 

512.5 

512.5 

512.5 

512.5 

512.5 

10.0 

535.5 

533.1 

530.9 

528. C 

525.6 

520.7 

516.5 

511.8 

507.0 

502*1 

597.3 

592.6 

588.3 

585.5 

581.1 

578.5 

576.6 

573.6 

575.7 

20.0 

551.C 

537.7 

533.1 

527.t 

521.2 

515. C 

506.3 

598.0 

589.3 

480*2 

571.2 

562.2 

553.5 

555.3 

538.3 

532.5 

528.1 

525.6 

525.6 

30*0 

531.0 

526.3 

519.2 

511.1 

502.5 

593.  C 

582.8 

572.0 

560.6 

448*4 

535.8 

522.9 

509.9 

397.3 

386.0 

376.1 

368.7 

365.3 

363.9 

40«0 

505.9 

598.il 

589.5 

579.3 

569.1 

558.5 

557.6 

535.5 

523.2 

409*2 

393.9 

377.5 

>60.5 

353.3 

327.1 

312.5 

300.5 

293.5 

292.6 

SOeO 

563.5 

556.2 

555.1 

533.1 

522.6 

512.9 

503.2 

392.5 

375.7 

364*9 

357.8 

328.6 

307.9 

286.3 

265.8 

255.7 

227.5 

216.6 

215.5 

60.0 

507.9 

399.7 

387.2 

375.2 

365.7 

359.2 

352.6 

355.0 

332.3 

317*8 

299.5 

278.3 

255.6 

229.2 

202.9 

177.. 

153.6 

137.6 

133.0 

70.0 

339.9 

331.1 

317.5 

305.3 

302.5 

301.5 

299.0 

293.i3 

283.5 

269*6 

251.1 

229.0 

203.5 

175.2 

155.0 

115.5 

85.6 

62.6 

55.6 

80.0 

261.7 

252.5 

238.1 

232.1 

238.3 

253.6 

255.5 

252.7 

235.0 

222*1 

205.5 

182.3 

156.2 

126.9 

95.5 

62.5 

29.3 

5.6 

0.0 

GG.O 

173.5 

166.0 

152.6 

165.1 

179.5 

185.5 

195.7 

195.3 

188.6 

177.1 

160.8 

139.8 

115.9 

87.0 

57.2 

26.8 

0.8 

C.O 

0.0 

A  and  B  are  the  azimuth  and  gradient  of  slope 
respectively;  in  the  programme  illustrated 
the  azimuth  is  given  from  0-  (north)  to  180* 
(south)  in  10  degree  intervals,  and  the  slope 
from  0*  to  90®  for  every  10®  degree  interval. 

W  is  the  hour  angle  e;^ressed  in  radians  ,  ,rith 
solar  noon  equal  to  0.000  radian;  negative 
values  are  before  noon,  and  positive  values 
Eure  after  noon;  in  the  programme  illustrated 
calculation  starts  at  0400  hours  cuid  terminates 
at  2000  hours;  the  calculation  is  made  at  20 
minute  intervals  (20  min.  =  0.873  radian). 


The  remaining  elements  of  the  programme  constitute 
calculations  from  the  formulae  used  to  compute  the  surface  values 
of  direct  solar  radiation.  These  elements  eure  found  in  the 
calculation  of  X(l),  X(2),  X(3) ,  T1  and  T2.  The  shadow  effects 
taken  into  account  are  those  caused  by  the  horizon  and  the  plane  of 
the  slope  itself.  In  the  programme,  QD  stands  for  the  cosine 
of  the  sun's  zenith  angle.  If  this  is  negative,  the  sun  must  be 
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DECLINATkON  OF  THE  SUN 
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Fig.  7  Seasonal  Oianoe  of  the  Latitude  of  Declination 

of  the  Sun 


TABLE  TWO 


Solar  Declinations  Selected  for  Cc»nputer  Programming 

I 


Period 

Declination 

r 

\xy/inxn; 

'.20  -  Jan. 23 

•  23*27' 

2.07 

Feb. 18;  Oct. 24  -  Nov.  29 

-15'’38' 

2.06 

Mar. 10;  Oct. 33  -  Oct. 23 

-07*49' 

2.04 

Mar. 31;  Sep. 13  *  Oct.  2 

00*00' 

2.02 

Apr. 21;  Aug. 23  -  Sep. 12 

07*49' 

1.99 

May  18;  Jul.26  ■*  Aug. 22 

15*38' 

1.97 

'  19  -  Ju.  25 

23*27' 

1.94 

If 
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below  tiie  horizon.  QT  denotes  the  cosine  of  the  angle  between  the  normal 
of  the  slope  and  the  direction  of  the  sun.  QT  will  be  negative  when 
the  sun  is  behind  the  slope.  Thus,  by  eliminating  QD  when  negative 
and  QT  when  negative  from  the  accumulation  of  the  values  ccmputed 
every  twenty  minutes,  it  becomes  unnecessary  to  determine  the  hours 
when  the  sun  shines  for  the  first  time  (t^  in  equation  (7)),  and  the 
last  time  (t^  in  equation  (7))  on  each  slope. 

It  is  apparent  that  the  computer  programme  illustrated  applies 
to  a  particular  latitude  (F)  and  refers  also  to  a  given  time  period  for 
which  the  angle  of  declination  of  the  sun  (FO)  and  the  extra-terrestrial 
radiation  (SCL)  remain  consteuit.  The  value  of  the  latitude  remains  a 
constant  which  fixes  the  location  of  the  area  being  excuoined.  The  values 
of  FO  2uid  SOL,  however,  will  veury  during  the  year. 

To  take  account  of  the  changes  in  FO  and  SOL  for  each  day  or 
week  of  the  yeeu:  is  time  consuming  and,  in  practice,  unnecessary  since 
the  influence  of  the  day  to  day  changes  in  the  sun’s  position  in 
equation  (1)  is  insignificant  as  conpared  with  the  influence  of  the 
other  terms.  Thus,  it  is  convenient  to  divide  the  year  into  periods 
for  which  a  siu table  declination  and  value  of  extra-terrestrial  radiation 
can  be  used.  Since  the  curve  of  the  latitude  of  the  svm's  declination 
against  time  is  sine-like  ,  the  procedure  followed  has  been  to  divide 
the  year  into  periods  of  unequal  length  for  which  a  given  declination 
and  value  of  can  be  regarded  as  a  true  mid-value.  Such  a  division 
is  given  in  Fig.  7  and  Table  Two.  The  length  of  the  periods  have  been 
chosen  bearing  in  mind  the  earlier  stipulation  that  an  accuracy  of  5% 

IS  appropriate  in  view  of  instrumental  limitat..ons.  A  division  of  the 
year  into  twelve  periods  has,  accordingly,  been  made , centering  upon 


Fig.  8  Regions  of  Barbados  for  Relief  San^ling 

X:  Couplicated  relief;  Y:  Intermediate  relief; 
Z:  Gentle  relief 


Gradient  Angle 


0  4  9  17  93  57  171 


Number  of  Sampling  pointe 


Fig.  9  Sampling  Points  and  Gradient  Angles  for  the 
Radiation  Happing  Grid  of  Barbados 

NOTE:  Bracketed  numbers  refer  to  the  sampling  areas 

in  Fig.  8 
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the  representative  declination  and  values  of  given  in  Tii>le  2. 

Mapping  the  Data 

In  order  to  laap  the  topograi^c  variations  in  stiort  wavs 
radiation  over  a  given  area  an  appropriate  base  nap  anst  ioe  pon^axed. 
I^is  involves  covering  the  area  of  study  with  a  grid  at  each  inter- 
secticm  of  which  is  iD2u:ked  the  aziauth  and  gradient  of  slope.  Ibe 
closeness  of  the  grid  must  be  dec?ded  in  terms  of  the  detail  of  the 
study  being  made  and  also  the  intensity  of  the  relief  of  tiie  area. 

Ihe  grid  system  used  for  the  map  of  Barbadoe  appearing  in 
this  report  was  made  fr<Mi  the  eighteen  sheets  of  the  1:10,000  topo¬ 
graphic  series,  prepared  in  1953  by  the  British  Colonial  Surveys 
(Basnayake,  1968).  these  map  sheets  have  on  them  a  basic  grid,  10  ca 
squcure.  Ihis  grid  provided  a  sufficiently  close  network  for  the 
southern  half  of  the  island,  with  its  gently-rolling  physiograpliy. 
the  more  hilly  central  and  eastern  parts  of  Barbados,  however, required 
a  closer  grid  system. 

To  determine  this  detail,  three  10  cm  squares  of  varying 
types  of  conplicated  relief  were  selected  (Fi% .  8) .  the  average 
gradient  for  each  square  was  calculated  using  successively  doubled 
numbers  of  sampling  points.  The  results  are  shown  ir.  Fig.  9.  the 
gradient  averages  reached  stable  vedues  at  between  17  and  33  points  per 
10  cm  square.  There fc^e,  the  figure  of  25  points  per  10  cm  square  was 
chosen  for  the  areas  of  most  complicated  relief  of  the  island,  this 
resulted  in  a  grid  pattern  of  2.5  cm  squares. 


An  identical  procedure  was  adopted  as  a  means  of  selecting 
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tbe  ^id  density  for  use  in  areas  of  intexaediate  relief.  As  a  result 
a  grid  of  5  cei  squares  uas  adopted  for  this  part  of  the  island.  Thus, 
it  was  possible  to  divide  the  island  into  the  three  regions  shown  in 
Fig.  8. 

TO  aeasure  the  gradient  of  the  slqpe  t  each  point  of  the 
resulting  grid,  the  nuiher  of  successive  contour  lines  over  standard 
distance  was  used  rather  than  the  difference  between  individual  a(?jacent 
contours.  This  procedure  was  followed  for  two  reasons:  it  is  easier  to 
count  the  successive  contours  ever  a  cosf/aratively  Icurge  pre-determined 
standard  distance  than  to  neasure  accurately  small  distances  between 
two  oontonrs;  tmd  to  use  a  staedaurd  distancp  for  greidient  measurement 
envies  an  average  value  to  be  obtained  that  takes  into  account  the  area 
arond  a  saspling  point,  rather  than  using  the  value  at  a  given  single 
point.  A  circle,  15  am  in  diameter,  was  used  for  gradient  analysis.  This 
proved  to  be  an  optimum  size  %diich  avoided  both  the  ambiguity  of  single 
contour  interpretation  arising  from  the  use  of  a  smaller  circle,  9nd. 
the  difficulty  of  determining  average  slope  direction  emd  counting  too 
maiy  contours  if  a  larger  circle  were  used. 

The  1:10,000  maps  used  to  prepare  the  base  map  have  contours  on 
them  at  20  ft.  intervals.  For  such  a  map,  therefore,  the  number  of  contours 
in  a  15  mm  distance  corresponding  to  different  angles  of  slope  is  given  by. 


15  tanX 
0.610 


(11> 


there  X  is  the  given  angle  of  slope  and  Y  the  number  of  contour 
crossings  in  15  mm  of  map  distance. 

The  calculation  of  the  values  of  azimuth  and  gradient  to  be 
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inserted  at  each  intersection  of  the  grid  was  effected  by  means  of  a 
clear  pletstic  protractor,  with  a  movable  pointer  of  semi-rigid  clear 
plastic  pivoted  at  the  origin  of  the  protractor.  The  pointer  was  marked 
with  a  central  longitudinal  line,  with  another  line  at  ri^t-angles 
to  it  and  passing  through  the  pivoted  origin.  A  circle  of  15  non  diameter, 
with  Its  centre  at  the  pivoted  origin,  was  also  marked  on  the  pointer. 

The  protractor,  with  the  pointer  attached,  was  placed  with  its  centre  of 
origin  on  the  selected  grid  intersection  of  the  topographic  map.  To 
obtain  the  eusimuth  of  the  s^ope,  the  pointer  vas  moved  So  that  its  mcuicer 
line  Weis  pexpendiculeir  to  the  major  contour  trend,  and  the  appropriate 
angle  was  then  read  on  the  circulair  margin  of  the  protractor.  For 
gradients,  it  was  simply  necessary  to  count  the  successive  contours  within 
the  15  mm  circle  and  to  read  the  corresponding  gradient  from  the  values 
in  a  table  prepared  from  equation  (11) . 

When  measuring  the  gradients  care  was  taken  to  ensure  that  the 
contour  lines  were  truly  successive,  and  were  not  repeating  the  same 
values  cis,  for  exanple,  in  the  case  of  those  doubling  back  on  themselves. 

In  the  case  of  open  valleys  or  spurs,  slope  measurement  was  done  on  the 
main  axis.  Narrow,  deep.,  valleys,  a  large  number  of  which  exist  in  the 
limestone  topography  of  Barbados,  were  ignored  in  the  slope  analysis. 

This  was  because,  in  general,  they  play  an  insignificant  role  in  the  regional 
picture  of  surface  radiation  variation,  and  especially  because  many  of 
them  were  found,  on  inspection  in  the  field,  to  be  choked  with  bushes  and 
trees  so  that  the  effective  surface  was  often  level  with  the  adjacent 
surroundings,  if  not  actually  higher. 

To  effect  the  irapping  of  global  radiation  by  using  the  base 
map  prepared  in  this  way,  a  computer  programme  was  made  which  calculated 
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the  value  of  the  gldaal  radiation  for  each  intersection  of  the  beise- 
map  grid.  This  calculation  is  made  frcat  the  radiation  data  of  a  single, 
representative  site,  «Aiich,  in  the  case  of  the  maps  appearing  in  this 
report,  constitute  the  observati<xis  of  global  and  sky-diffuse  radiation 
made  in  1964  at  the  Brace  Experiment  station,  in  the  grounds  of  the 
Bellairs  Research  Institute,  St.  James  (Tout,  1967,  1968). 

For  mapping  purposes  two  sets  of  data  Ccurds  were  prepared: 
cme  showing  the  azimuth  and  the  other  showing  the  gradient  for  each 
intersection  of  the  base-mc^  grid.  Ihe  maximum  possible  number  of 
columns  per  card  for  this  purpose  was  twenty,  since  four  units  cure 
needed  to  indicate  the  azimuth  in  degrees  and  whether  it  is  positive  or 
negative.  Thus,  the  grid  covering  Barbados  Weis  grot^d  suitably  into 
major  divisions,  each  division  being  twenty  grid-  points  wide.  The 
progranme  thus  calculated  the  global  radiation  values  at  the  grid  points 
in  one  continuous  length,  twenty  columns  wide.  The  resulting  print-out 
was  then  divided  suitably  so  as  to  conform  to  its  actual  pattern  on  the 
map.  To  ensure  the  correlation  of  points  on  the  map,  an  appropriate 
overlap  was  made  in  the  division  of  the  grid  into  twenty-column-wide  units. 

The  calculations  effected  by  the  computer  are  illustrated  in 
Fig.  10.  The  latitude  (F)  is  expressed  in  radians.  After  reading  the 
cards  giving  the  azimuth  and  gradient  of  each  point  in  degrees,  the 
values  at  these  points  are  converted  into  radians  in  steps  12  and  13. 

Step  15  reads  the  year,  month,  and  day  (LI),  the  declination  of  the  sun 


In  practice,  it  Weis  found  convenient  to  confine  the  division  into  twenty 
columns  to  the  region  of  complicated  relief  on  the  island.  The  remainder 
of  the  island,  where  the  base  grid  was  wider,  was  divided  into  fourteen 
column-wide  units,  and  the  computer  programming  for  calculating  this  part 
of  thi  island  was  modified  accordingly. 
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C  GLOBAL  RACUTICN  MAP  BAReACCS 

CIPEK$ICK  All20i20l •8(120*201. XI3>>SI 120,201 
C  INPUT  LATITUDE  IN  RAOiAhS 
F-C.22C4 

C  READ  THE  AUMBER  OF  CARDS 
REA0(5,88ei  L 
888  F0PNATII3) 

C  READ  AZIHUTH  ANC  GRACIEAT  CF  SLOPES 

REAO(5*100I  i(A<I*JI,  J«l,20t,I>l*L> 

100  FORMAT  (20F4.0t 

REACI5,101i  (tedfJt  •J<1*2C>*  1>1*LI 

101  FORMAT  (2C-FS.0) 

CO  880  i*l,L 

00  880  J>1.2C 

IFIBtt tJI.EC.4C.CI  6C  TC  €80 
A(l.Ji«A(I,Ji*0.Cl7<S 
B(I,J)«B(I,J)»0*CI745 
880  CONTINLE 

661  REAOI5.1C2,EKO»/44»L1,  rC,  IF,  CIR,  OIF 

102  FORPAfde,  F4.0.  12,  2F6.Ct 
liRITE  (6.456) 

456  FORMAT  IlHlI 

hRITE  (6,103IL1,  EC.  IF,  CIR,  CIF 

103  F0RMAT(6H  DATE#  •  18,  JSf"  CECLlNATtON*  ,F4,C,  I2//14H  DIRECT  INSOL 
lATIONA  ,F6, l/24h  SKV-DIFFLSE  IKSClATICNA  ,F5.1) 

FIF-IF 

FIF=FIF*0,0167 
IF(FO)  2C,  20,  21 

20  FC»(FC-FIF)*0.01745 
GC  TC  22 

21  F'-*(FCaFIF)*0.01745 

22  G*0.50 
134  SS-0.0 

tix-2,16 

155  CO«COS(FCI*CCS(r)*COS(fc)*SIN(FO)PS  MF) 

G1*1.0/ICCaO.OOCI1) 

IF(AeSICl!.GT.4(i.0)  C1>3C.C 
|FIQ0,LT.0,0)  C0>0,0 
ss«ssAeo,oAGA*ci*co 
l,«N+0,1746 

IF(H,LT,2.Z!  go  TC  155 
C  COMPARISON  CETmEEN  CIR  ARC  SS 
0-CIR-SS 

IFIABSICI-  5.01  162,  162,  161 

161  IFIDIR.GT.SSI  GO  TO  163 
G»G-0,C5 

GO  TC  154 
163  G-G+0.C15 
GO  TO  154 

162  RRITE  16,  1711  G 

171  FORMAT  (17H  TRANSMISSIVITY#  , 

00  778  I»1,L 
CO  778  J=1,2C 
IF(B(I,JI.EC,90.)  GO  TO 


Fig.  10  The  Computer  Programme  for  Calculating 
Global  Radiation  over  Barbados 


(ccsitlnued  over  page) 
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0049 

xiii»cos{Aii«jn*si6ieiif  j») 

0090 

XI  2I>S  IKIBI 1  >31  }6S1NU(  I  ,Ji  1 

0091 

y(3l>CCS(BCItJlt 

0092 

Tl«(Xll|6SINIFI-»X(3i*CCS|FM*CCSIF0l 

0053 

T2-t-Xll|6C0SIF)«X(3i*SI6IF) 1*5101 FCi 

0094 

6»-2«16 

0055 

Slit Jl*0.0 

0054 

5 

Q0*C0SCFC|6G0S(F»«CCSIbMSIMlFC>*SlNIF) 

0097 

Ci«1.0/{C0*0.0001) 

0090 

IFI01>2.ei  121>  121>  122 

0059 

122 

IFIQl.eE.2.91  01«2.90 

0060 

IF(01.eE.3*0>  01*3.05 

0061 

IFI01.eE.3.t)  Cl-3.21 

0062 

1F(01.€E.3.3I  01*3.39 

0063 

lFi01.CE.3.SI  01*3.59 

0064 

IFIOl.EF.3.7)  01*3.82 

0065 

iF(01.CE.4.0t  01-4.07 

0066 

IFIOI.CE.4.61  01*4.72 

0067 

1F(Q1.£E.5.0I  01*5.12 

0060 

IFIOl.CE.5.51  01*5.60 

0069 

IFIQI.CE.6.0)  01-6.18 

0070 

IFIQl.CE.  6.51  01-6.68 

0071 

IFIOl.EE.  8.C»  01*7.77 

0072 

IF<C..£E.  8.81  01*8.90 

0073 

IFIOl.OE.lO.CI  01*10.39 

0074 

IFIOl.OE. 13.01  01*12.44 

0079 

IFIOl.EE. 17.01  01*15.36 

0076 

!FI01.CE. 24.01  01*19.79 

0077 

121 

IF(A8S(ail.GT.45.01  01*30.0 

0070 

02*-X(21«S!NfU16CCSfF0t6Tl«CCSItfl 

0079 

OT<«a2*T2 

0060 

I F 10  C. IT.O.O.OX.OT.LT. 0. 0 1  Ot*C. C 

0081 

SI  I > J 1 *5 U t J  1680.06C6*C1*GT 

0082 

N*^0.  1745 

0083 

iFlh.LT.2.21  GC  TO  5 

0084 

SII>Jt«SII>Jl«0iF*{CCSI8II.JI/2.Cl >•*2 

0085 

eo  TO  778 

0086 

414 

SI  I.JI>0.0 

0087 

778 

CONTIMF 

0088 

WRITE  {U.6321  IISiI><Of  J*l>201.I*l>LI 

0089 

632 

F0R867I2CF5.01 

C090 

GO  TO  661 

0091 

999 

STCR 

0092 

END 

Pig.  10  (Continued) 
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(FO,.  IF),  and  the  direct  (DIR)  and  sky-diffuse  (DIF)  radiation  for  the 
day  in  question.  The  declination  of  the  sun  is  e3q>ressed  in  degrees  (FO) 
and  minutes  (IF) ,  and  this  value  is  converted  to  r€idians  in  steps  21  to 
2C.  The  direct  and  sky-diffuse  radiation  are  daily  totals  in  langleys. 
Thus  the  programme  can  be  used  for  a  particular  day,  or  else  for  a  longer 
f ')riod  by  using  the  mean  daily  value  for  the  period  in  question.  Steps 
26  to  43  effect  the  calculation  of  atmospheric  transmissivity  (G). 

This  is  achieved  by  first  calculating  a  daily  total  (SS)  on  the  bcusis  of 
a  transmissivity  of  0.50  and  accumulations  every  40  minutes.  This  total 
IS  then  compared  with  the  direct  radiation  total  alrecufy  sillied,  and  a 


re-calculation  eiltering  the 


'f  G,  either  by  decreasing  it  (step 


40)  or  increasing  it  (step  42. 


fected  until  the  Vcdue  of  SS  is 


within  5  langleys  of  the  value  of  measured  direct  radiation.  The  value  of 
G  which  achieves  this  is  then  used  to  calculate  the  direct  radiation  at 
each  grid  point  in  a  manner  similar  to  that  already  ej^lained  in  connection 
with  the  base  programme.  The  contribution  of  sky-diffuse  radiation  is 
added  in  step  84. 

The  programme  illustrated  in  Fig.  1C  is  arranged  for  a  print¬ 
out  in  table  form  (Fig.  11) .  This  print-out  gives  the  values  of  globed, 
radiation  at  each  grid  intersection  of  the  base  map.  Plotting  the  data 
eind  drawing  the  relevcint  iseurithms  by  hand  thus  becomes  routine  carto¬ 
graphic  procedure.  Although  this  perhaps  sounds  an  unnecessarily 
lengthy  way  to  produce  the  final  map  in  view  of  the  availability  of 
computer  mapping  programmes,  in  fact  it  was  found  the  most  appropriate 
way  for  making  the  maps  of  Barbados.  This  conclusion  was  reached  after 


trying  out  the  computer  mapping  programme  illustrated  in  Pig.  12. 
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Fig.  11  A  Sanple  of  Gldt>al  Radiation  Values, 
Barbados 
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C  (t«CIATION  STUCVt  GLOBAL  RACIATION  MAP*  BARBACOS. 
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Fig.  12  Con«>uter  Programme  to  Map  Global 
Badiation;  Barbados 
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The  print'out  for  this  progranone  (Fig.  13)  uses  a  syndbolism 
based  on  an  isairithmic  interval  of  2S  langleys.  Tlie  arrangeaent  of 
this  symbolism  is  based  on  groups  of  three  letters  separated  by  blaztks. 
Since  dedly  totals  of  global  radiation  cure  unlikely  to  be  below  100 
langleys  in  the  tropics,  the  first  symbol  (*)  indicates  a  range  of 
75-100  langleys.  Hie  other  values  then  follow  by  letters  of  the  alphabet: 

*  ABC  DBF  GHI  JKL  MNO  PGR  STU 

From  this  sequence  it  follows  that  A  refers  to  125-149  langleys,  B  to 
150-174  langleys,  C  to  175-199  langleys,  the  blank  between  C  and  D  to 
200-224  langleys,  cnl  so  on<.  Hie  plotting  scale  on  the  computer  map  is 
half-an-inch  between  each  grid  point,  and  the  computer  effects  an  inter¬ 
polation  between  the  calculated  values  at  each  of  these  points. 

A  sample  of  a  map  of  pcurt  of  Baurbados  prepared  in  this  way, 
with  isarithms  drawn  in  by  hand,  is  given  in  Fig.  13.  The  map  illustrates 
the  advantages  of  the  system  as  a  rapid  means  of  providing  a  general 
picture  of  the  distribution.  For  a  map  of  the  whole  of  Barbados,  however, 
modification  of  the  programme  is  necessary  in  view  of  the  different 
density  of  grid  points  used  to  sanple  the  varied  topography  of  the  island, 
or  alternatively  a  rather  conplicated  arrangement  of  enlargement  euid 
reduction  of  the  maps  of  different  regions  has  to  be  done  to  bring  them 
together  on  one,  suitable  scale.  Even  with  such  modifications,  all  of 
which  were  tested,  the  resulting  map  presented  a  rather  generalized 
picture.  After  a  number  of  trials  of  different  systems,  therefore,  it 
was  decided  to  use  the  digital  print-out  illustrated  in  Fig.  10  as  the 
most  appropriate  method  for  preparing  the  maps  illustrated  in  this  report. 
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Hourly  Radiation  Mapping 

AlthoL.^  the  aethod  for  evaluating  ttie  surface  vaxiatioBS  of 
glcdaal  radiation  which  has  been  outlined  is  intended  priaarlly  for  use 
with  daily  totals  or  longer  periods,  the  tests  for  accuracy  giusn  to  the 
system  at  Mont  St.  Hilaire,  Quebec,  suggest  that  it  can,  in  fact,  be 
used  for  shorter  periods,  with  prdoably  an  hour  as  the  miniawi  reliable 
time  period.  It  is  useful  in  the  tropics  to  examine  the  radiation 
distribution  over  such  a  short  period.  This  is  because  the  sm  rises 
quickly  and  its  heating  effects  are  strong  sotm  after  sunrise  and  until 
just  before  sunset.  By  the  time  the  sun  has  reached  an  altitude  of  30° 
above  the  horiz<xi  topographic  inflvences  becoee  scaeidiat  reduced.  The 
varied  pattern  of  heating  in  the  early  morning  hours,  however,  can 
theoretically  influence  the  develcpment  of  early  morning  convection. 

In  Barbados  the  steepest  slopes  on  the  island  teitd  to  be  generally  towards 
the  east.  Thus,  xt  was  felt  desirable  to  arrange  a  conputer  programme 
to  calculate  the  topographic  Vciriations  of  global  radiation  on  an 
hourly  basis. 

Such  a  programme  is  illustrated  in  Fig.  14.  It  is  basically 
the  same  as  those  shown  in  Figs.  10  cutd  12,  but  differs  sometdiat  in  the 
data  supplied  tc  it.  These  data  are  indicated  in  step  15  of  Fig.  14. 

They  consist  of  the  year,  month,  and  day  (LI),  the  angle  of  declination 
of  the  sun  (FO,  IF) ,  tha  extra-terrestrial  radiation  (SOL)  in  langleys 
per  minute,  the  atmospheric  transmissivity  (G)  for  the  hour  in  question, 
the  sky-diffuse  radiation  in  langleys  for  the  hoxir  (DIF) ,  and  the  time, 
which  is  given  as  eui  hour  angle  in  degrees  and  decimals  of  a  degree, 
with  solar  noon  as  zero  and  negative  values  before  noon.  This  time  vadue 


1 

■  ? 

-  40  - 

a 

C  PA'-IAT  rr  STUfY  CF  FAPEfOCS  HCOFLV  M4PPIM0 

• 

0001 

01PLK5 IC4  «< 120t2CI .61 12C.20l*X<3t.SI 120.201 

1 

C  IfJPLT  LATITuOc  in  PACUAS 

• 

0002 

F*C.22t9 

‘ 

C  RF*C  TM6  60PP?P.  OF  C6RCS 

0003 

RF«c(5.d<i8l  L 

0004 

888  FCFMXICUI 

C  RFC.  C2IPUTH  CNC  GR.>CICAT  CF  SLOPFS 

0005 

PFC0(5.100I  ((CII.JI.  J<1.20t.I>l.L> 

J 

C0C6 

100  fi-pnCT  l2CF4.0i 

i 

0007 

RFCKS.lOli  t(6(I,JI.J«1.2CI.  I>1*LI 

oocs 

101  FCRPCT  (20F3.Ot 

coco 

CO  880  I«l.t 

0010 

no  880  Jp1.20 

0011 

IFl4(I,Jl.eC.9C.0l  GC  TO  660 

0012 

*(  ItJ|pC(I.JI»'0.C1745 

0013 

PlI.JMBd.JlPO.ClTAS 

0014 

880  CONTIACF 

[ 

i 

C015 

661  RFCp(5.1C2t640«599ILl>  FC»  IF.  SCLt  6.  DIF,  TIPF 

1 

0016 

102  FORHCTlie.  F4.C.  12,  4F6*C> 

1 

i 

0017 

k-TMCPO.OlTAh 

0018 

WRire  (6,9561 

0019 

956  F0(-8AT  (IHlt 

1 

0020 

hRIT6(6.103l  Ll.FC.  IF,  C,  DIF,  TINr 

1 

t. 

0021 

103  F0RPAT(6H  0AT6R  .  18,  136  CECilAATICN*  ,F4.0,  I2//1TH  TRANSPISSIVI 

j 

ITVA  ,F6«3/24H  SKV>D1FFI.SE  INSGLATlCNf  ,F5.1/6H  TIME*  ,F6.2I 

i 

i  ^ 

0022 

FIF-IF 

1 

0023 

FIF-FIF»C«0167 

0024 

IFlFCI  2C,  20,  21 

■1 

0025 

20  FC-IFC-FIFi*C.01745 

H 

0026 

GO  TC  162 

0027 

21  FC«(FC4FIF|*C.C1745 

0028 

162  CONTIKLE 

0029 

DC  776  I-1,L 

0030 

CO  778  J«1.20 

0031 

IF(8(I,JI.EC.90.t  GO  TC  414 

0032 

X(li— rCSIA(I.J|»*SIA(6(!,jn 

0033 

X(2l>SIMetI,Jil«SIN(A(I,J»l 

0034 

XI3)>CCSIB(I,JII 

• 

i 

• 

0035 

T1«(X(  1I4SIMF  MX(3I*CCS(FI)*CCSIFCI 

- 

1 

0034 

T2»(-X(1»*CCS«F:4XI31*SIMF}16S1MFC» 

i 

0037 

S(I,JI>C.C 

1 

i  i 

0038 

5  QC>CCS(FCI4CCS(F(*CCSIbl«SlN(FOI*SIN|FI 

■'  i  ! 

0039 

01>1.0/(G040«0001I 

‘ 

0040 

IF(01>2«8)  121,  121,  122 

1  i 

0041 

122  1F(C1.EE.2,9»  C1-2.9C 

i 

0042 

IF(01*CE.3.CI  GIO.OS 

i  1 

0043 

IF(Cl.r-E.3.1l  Gl>3.21 

C044 

IF(Q1.EE.3.3I  Gl«3.39 

*  . 

1 

0045 

IF(01.CE.3.SI  C1-3.S5 

0046 

IF(01.CE.3.7I  Cl«3.82 

0047 

IFtOl.CE.4.01  Cl«4.07 

, 

0048 

IF(01.GF.4.6I  Cl-4.72 

0049 

IFIOl.Gc.S.C)  Cl«5.12 

5  J 

0050 

IF(0'l.CF.5.5l  C1>6.6C 

0051 

IFIOI.Ct.e.OI  Cl*6.18 

1 

0052 

IF(01.Gc«  6.51  Cl*6.88 

0053 

IFIOl.GF.  6.0  Cl-7.77 

0054 

IFICl.GE.  8.81  Cl«8.90 

0055 

IFfOl.CE.lJ.OI  C1-1C.35 

0056 

tF(Ql.GE.13.CI  Cl-12.44 

. 

0057 

IFIOl.GF.ir.CI  Gl>15.36 

5 

0058 

IFIOI.CF.24.01  Cl-19.75 

• 

. 

0059 

121  IFIABSICII. 87.45.01  C1-30.C 

:  : 

0060 

02>-X(2|4SI«l(i,)*CCS(FCl4ll«CCS(k» 

0061 

CT»g2*T2 

0062 

IF(0C.IT.C.C.CR.C7.LT.C.C»  0T>C.C 

0063 

S(  I .  Jl  -S 1 1 ,  J  U60.0«SCL«6*«C1«0T 

0064 

S(I.JI«!,'I.Jl4PlF*(CC$(e(I.JI/2.CI»**2 

0065 

GO  TO  776 

C06i 

414  SI  I. Jl <0.0 

0067 

778  C06TIMF 

- 

0068 

kRITE  (6.6321  ((Sd.jlt  J>t .201  .I>1.LI 

0069 

632  FORMAT (2CF5.0I 

0070 

CF  TC  661 

0071 

999  STCP 

0072 

E»r 

. 

Fig.  14  Campatex  Pzogrmmt  far  Calcolatiiig  accrl/ 

f 

1 

Global  Kadiatios  in  Barbados 

I 

1 

_ 

A 

-41 


is  then  converted  to  radians  in  step  17. 

Ihe  value  of  sky-diffuse  radiation  is  supplied  directly  for 
the  hour  in  question  frcnn  the  radiation  observations  of  the  sita  being 
used.  To  obtain  the  atmospheric  transmissivity  (6)  it  is  necessary 
first  to  obtcdn  the  height  of  the  sun  above  the  horizon  for  the  given 
hour.  Uiis  is  provided  by  a  programme  of  vdiich  a  print-out  of  resiilts 
is  illustrated  in  Fig.  15.  In  this  exauple,  the  solar  declination  (FO) 
is  for  June  21st  and  shows  that  at  6  am  solar  time,  hour  angle  (W)  =  0.00, 
the  hei^t  of  the  sun  is  5.04°  above  tl.'e  horizon  in  Barbados.  Hie 
values  for  W  are  given  every  twenty  minutes,  from  which  the  height  of 
the  sun  at  the  mid-point  of  the  relevant  hour  can  be  readily  found. 

The  results  of  a  second  programme  (Fig.  16)  are  then  used  to  obtain  G 
by  reference  to  the  height  of  the  sun  above  the  horizon  (EJ)  and  the 
value  of  the  direct  radiation,  in  langleys  per  minute  (Q) ,  observed  at 
the  station  during  the  ht,  ur.  llTiese  data  are  then  inserted  into  the 
computer  programme  and  the  calculation  of  the  global  radiation  at 
each  grid-point  of  the  base  map  proceeds  in  the  usual  way,  with  a 
digital  print-out  for  results  (Fig.  17) . 

It  will  have  been  noticed  that  neither  in  the  mapping  of 
hourly  values  nor  of  daily  totals  has  any  allowance  been  made  for  the 
influence  of  local  sky-line  c^struction  on  the  topographic  variation 
of  radiation.  A  method  to  allow  for  this  %dien  necesscury  has  been 
devised  cuid  incorporated  into  the  computer  calculations  (CXunira,  1969) . 
However,  an2LLysxs  shows  that  the  influence  of  a  local  sky-line  is 
negligible  on  daily  totals  except  in  mountainous  areas  of  middle  and 
hi^.  latit*udes,  and  for  hourly  values  only  when  the  sun  is  low  and  the 
local  relief  particularly  rugged.  It  needs  to  be  allowed  for  in  the 
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tropics,  therefore,  only  under  exceptional  circumstances,  %dien 
for  exanple,  a  detadled  study  over  short  time  periods  is  underway. 

In  Beurbados,  the  local  relief  is  nc%diere  so  strong  as  to  influence 
materially  the  distribution  pattern  of  global  raidiation  at  the  scale  of 
the  present  study.  Consequently,  no  allowcmce  has  been  made  for  such  an 
influence  in  the  work  discussed  ::n  this  report. 


equipKDt  at  terace  Experiaental 
Statioo.  Upper  photo:  Kipp  solariaeter. 
Lower  photo:  .shielded  silicon  cell,  vith 
r  ia|kbel  1  'Stokes  sunshine  recorder  nearby. 
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IV 

SELECTED  GLOBAL  RADIATION  MAPS  OF  BARBADOS 


The  methods  of  calcu  .tion  and  the  computer  programming  tdilch 
have  been  described  have  been  used  to  nuike  a  number  of  global  radiation 
maps  of  Barbados.  These  cure  somewhat  preliminary  in  nature,  and  it  is 
not  the  purpose  of  this  report  to  offe"  either  a  detailed  analysis  of 
them  or  a  substantive  scientific  interpretation  of  the  results  which 
they  show.  Some  genered.  comments  on  a  few  of  the  maps,  however,  may 
serve  to  illustrate  further  the  scope  and  usefulness  of  the  procedures 
which  have  been  developed. 

All  the  data  for  the  maps  illustrated  here  were  derived  from 
observations  made  in  1964  at  the  Brace  Experiment  starion,  located 
in  the  grounds  of  the  Bellairs  Research  Institute  on  the  west  coast  of 
the  island  (Tout,  1968).  An  Eppley  pyranometer  and  a  Kipp  and  Zonen 
pyrcuiometer,  both  calibrated  at  the  Canadian  Meteorological  Heculquarters 
in  TOront'' ,  were  used  interchangeably  for  measuring  global  radiation. 

A  silicon  cell  unit,  with  a  shadow  band,  Mas  tised  for  the  skydiffuse 
component.  Corrections  were  made  to  allow  for  the  fraction  of  the 
radiation  screened  by  the  shadow  band  itself,  and  the  resulting  values 
were  subtracted  from  measured  global  radiation  values  to  obtain  the 
value  of  direct  solar  radiation. 

The  instrments  were  connected  to  Honeywell  strip  chart 
recording  potentiometers,  the  continuous  record  thus  obtained  being 
analysed  by  hand  zo  obtain  hourly  values.  There  was  scae  obstruction 


Pig.  18  Mean  Daily  DLstrix>ntioa  of  Global  Kadiation 
in  Bazbados,  1964 
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to  the  site  by  trees.  This  influenced  daily  totals  by  no  more  than  a 
maximum  of  3%  (Tout,  1967) .  However,  the  presence  of  the  obstruction 
means  that  the  rrecords  near  the  beginning  emd  end  of  each  day  are  unreliable 
for  hourly  analyses. 


TABLE  TliREE 

Mean  Daily  Radiation  Totals  for  Sele :ted  Months ,  1964 
(Ictngleys  per  day) 


Month 

Direct 

Sky-Diff’-ise 

Global 

Jidy 

241 

210 

451 

September 

212 

199 

411 

December 

318 

114 

432 

The  first  three  maps  used  as  illustration  (Figs.  18,  19,  20) , 
show  the  mean  d2d.ly  value  of  surface  variation  in  global  radiation  for  July, 
September,  and  December,  1964.  It  is  interesting  to  note  from  Table 
Ohree  that  the  mean  deiily  valve  of  globed,  radiation  recorded  for  each  of 
these  months  was  rather  similar.  This  similarity  fror  one  month  to 
another  is  a  feature  of  global  radiation  <^servations  in  the  tropics 
when  such  measurements  etre  confined  to  instruments  well-exposed  in  a 
horizontal  position.  Figs.  18,  19,  and  20,  however,  shew  that  such 
similarity  successfidly  hides  the  topographic  variation  %du.ch  actually 
exists.  This  variation  is  expressed  in  the  mc^s  illustrated  partly  by 
t:i3  range  o£  values  over  the  island,  and  partly  by  the  actual  spatial 
pattern.  The  former  are  low,  approximately  30  langleys  per  day,  din  the 
two  nigh  stm  months  of  July  and  Septexher,  but  the  range  exceedls  100 


1 


i 
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langleys  per  day  in  December.  There  is,  however,  more  spatial  variety 
in  the  July  map  than  in  the  September  map.  This  may  be  explained  by 
the  fact  that  the  sun  is  further  from  the  zenith  at  noon  over  Barbados 
in  July  than  in  September,  especially  during  vhe  early  part  of  the 
montlu  The  relatively  low  range,  however,  reflects  the  conparatively 
large  contribution  of  sky-diffuse  radiation  to  global  radiation  totals 
(Table  Three)  since  sky-diffuse  radiation  is  not  as  sensitive  as  direct 
solar  radiation  to  topographic  influence. 

In  December,  1964,  some  75%  of  the  global  radiation  total  was 
direct.  With  the  sun  at  its  lowest  elevation  of  the  yecir  during  this 
month,  the  resulting  spatial  pattern  expressed  considerable  diversity. 
The  influence  of.  surface  geometjry  on  the  spatial  distribution  of  the 
direct  radiation  component  under  these  conditions  is  accordingly 
manifested  in  the  large  range  of  daily  values  from  place  to  place  during 
this  month. 

A  specific  exanple  of  the  influence  of  solar  elevation,  even 
in  the  tropics,  when  the  direct  component  in  global  radiation  totals  is 
high,  is  illustrated  in  figs.  21  and  22.  The  measured  direct  radiation 
totals  on  the  two  days  shown  ware  almost  identical :  438  langleys  on 

June  18,  and  436  langleys  on  December  28.  Sky-diffuse  radiation  for  the 
two  dates  was  152  langleys  and  74  langleys  respectively.  The  resulting 
totals  of  global  radiation,  590  langleys  for  June  18  and  510  langleys 
for  December  28,  are  typical  of  the  general  range  of  radiation  values 
in  the  tropics  when  measurements  are  confined  to  the  horizontal.  A  much 
greater  contrast  is  apparent  viien  the  element  of  t<^>ographic  variation 
is  introduced. 

In  comparing  Fig.  21  with  Fig.  22  it  must  be  noted  that  the 
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scale  of  shading  differs  from  one  map  to  the  other.  Hie  June  map  has 
an  isarithmic  interval  of  25  langleys  and  that  for  December  of  75 
langleys.  The  cartographic  problem  of  showing  adequately  the  two 
patterns  by  the  same  isarithmic  interval  at  the  scale  permissible  at  the 
page  size  of  this  report,  itself  underlines  the  contrast  between  the 
surface  variations  of  global  radiation  on  the  two  dates.  The  range  of 
values  on  June  18  weis  little  more  than  75  langleys;  that  for  December 
28  exceeded  225  langleys.  Thus,  an  interval  of  75  langleys  for  the 
June  map  would  have  suggested  a  topographic  uniformity  which  did  not, 
in  fact,  exist,  while  a  25-  langley  interval  for  tlie  December  map  would 
have  produced  a  complexity  too  great  to  be  shown  on  a  small-scale  map. 

These  comments  alone  en^hasize  the  contrasts  in  the  topographic  variation 
of  global  radiation  during  the  two  days  —  two  days  with  generally  similar 
solar  radiation  totals,  and  an  atmospheric  tramsmissivity  and  pattern  of 
cloud  chcuracteristic  of  a  clecir  Barbados  day. 

The  remeiiring  n^aps  illustrate  mapping  on  the  basis  of  hourly 
radiation  values.  An  example  for  a  specific  date,  June  18,  1964,  is 
given  in  Figs.  23  and  24.  The  morning  hour  chosen,  0800-0900  hours,  is 
the  first  hour  of  the  day  in  which  the  observing  instruments  were  clearly 
unaffected  by  the  obstruction-^  referred  tc  earlier  eis  existing  at  the  site. 
The  midday  hour  was  chosen  for  con^jarison  becau'-e  the  sun  was  then  at 
Its  highest  elevation  of  the  day. 

The  day  in  question  was  clear  to  begin  with,  but  clouds  developed 
later.  Thus  of  the  47  lemgleys  recorded  on  a  horizontcd.  surface  between 
0800  and  0900  hours,  ,38  were  contributed  by  direct  radiation.  By  contrast, 
one-third  of  the  70  langleys  recorded  between  1200  and  1300  hours  were 
contributed  by  sky-diffuse  radiation.  Whereas  atmospheric  transmissivity 
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was  0.80  in  the  rooming  hour  it  had  been  reduced  to  0.63  by  midday. 

Con^eurison  of  the  two  maps  shows  clearly  tl.e  greater  spatial 
diversity  at  the  rooming  hour  cuid  also  the  higher  radiation  income 
of  the  east- facing  slopes  than  at  midday.  Indeed,  some  areas,  such  as 
northwest  of  Belleplane,  chemged  little  in  the  value  of  their  global 
radiation  Income  between  the  rooming  hour  and  midday,  whereas  the  rest 
of  th's  island  experienced  an  increcised  radiation  total  as  the  sun 
approached  its  highest  point. 

The  monthly  maps  of  mean  values  in  December  (Figs.  25  and  26) 
also  shows  a  tendency  for  east-facing  slopes  to  have  a  relatively  high 
global  radiation  income  in  the  morning  hour.  That  is  not  true,  however, 
of  the  two  July  maps  (Figs.  27  and  28).  The  difference  may  be  related 
to  the  contrasts  in  the  contribution  of  direct  radiation  to  global  hourly 
totals  in  the  two  months  (Table  Four) .  Whereas  in  December  thi  Tiean 
daily  contribution  of  direct  radiation  was  76%  at  0800-0900  hours  and 
73%  at  1200-1300  hours,  in  July  the  contributions  at  the  same  hours  were 
both  57%. 


TABLE  FOUR 


SELECTED  MEAN  HOURLY 

RADIATION  TOTALS, 

1964 

Month 

Time 

Direct 

SHy-Diffuse 

Global 

July 

0800-0900 

20 

15 

35 

1200-1300 

34 

26 

60 

December 

0800-0900 

26 

8 

34 

1230-1300 

44 

16 

60 

Fig.  24  Distribution  of  Global  Padiation  in  Barbados 
1200  -  1300  hours,  June  18,  1964 
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Fig.  27  Mecin  Value  of  Global  Radiation  Barbados 
0800-0900  hours,  July,  1964 
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V 

CONCLUSION 

The  aim  of  t]iis  report  has  been  to  outline  and  expletin  methods 
which  have  been  developed  to  portray  the  topographic  variaticsis  in 
global  radiation  by  using  the  radiation  cd>servations  of  a  single, 
representative  site.  The  sample  maps  used  as  illustration  provide  an 
indication  of  what  the  method  reveals  in  a  tropical  context.  Brief 
and  superficial  as  cire  the  comments  made  on  the  latter,  it  would  appear 
that  the  principal  conclusion  which  can  be  drawn  is  that,  even  in  the 
tropics  and  in  an  area  of  generally  ly  moderate  relief,  the  actual 
surface  variation  of  global  radiation  totals  and  intensities  are 
greater  than  the  normally  made  observations  imply.  These  variations 
not  only  provide  an  important  surface  pattern  at  a  given  time,  but  also 
vary  in  their  character  at  different  times  of  year  and  under  different 
weather  conditions.  In  cin  atmosphere  as  sensitive  to  disturbance  as 
the  humid  atmosphere  of  the  Barbados  region,  small  variations  in  the 
surface  radiation  may  well  have  important  co..sequences  for  convective 
development  through  the  day.  For  this  reason  alone  a  fuller  investigation 
of  the  pattern  of  short-wave  radiation  income  revealed  by  the  preliminary 
maps  of  this  report  seems  worthwhile. 

There  are, however,  wider  implications.  Many  aspects  of  economic 
development  are  related  to  a  knowledge  of  surface  radiation  balance.  This 
is  particularly  so  in  relation  to  problems  of  agricultural  development, 
plant  growth,  and  water  supply.  To  understand  better  the  topographic 
variations  of  the  surface  energy  budget  is  a  firsi-  step  towards  a  better 
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understanding  of  the  ecological  envixonaent.  the  diversity  of  this 
environaent  is  intiaately  related  to  the  diversity  revealed  by  this 
e^licaticm  of  the  analytical  and  naming  aethods  outlined  in  Ais 
report.  Tropical  envixonaents  are,  indeed,  not  nearly  sc  unifozB 
as  one  is  often  led  to  suppose. 
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